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' We investigate hard-thermal-loop (HTL) corrections to the CP-asymmetries in neutrino 

and, at high temperature, Higgs boson decays in leptogenesis. We pay special attention to 
the two leptonic quasiparticles that arise at non-zero temperature and find that there are four 
contributions to the CP-asymmetries, which correspond to the four combinations of the two 
leptonic quasiparticles in the loop and in the final states. In two additional cases, we approximate 
the full HTL-lepton propagator with a zero-temperature propagator that employs the thermal 
lepton mass me{T), or the asymptotic thermal lepton mass ^/2mi{T). We find that the CP- 
asymmetries in the one-mode approaches differ by up to one order of magnitude from the full 
(— I ' two- mode treatment in the interesting temperature regime T ^ Mi. The asymmetry in Higgs 

boson decays turns out to be two orders of magnitude larger than the asymmetry in neutrino 
decays in the zero-temperature treatment. The effect of HTL corrections on the final lepton 
asymmetry are investigated in paper H of this series. 



CN ■ 1 Introduction 

The matter-antimatter asymmetry of the universe is usually expressed as 

= (6.16 ±0.16) X 10"^°, (1) 



riB - 
r] = 

^ i which has been inferred from the 7-year WMAP cosmic microwave background (CMB) anisotropy 
5-( I data [1] , where ub, nB, and are the number densities of baryons, antibaryons, and photons, 
and the subscript implies present cosmic time. Leptogenesis [2] is a very attractive model in 
this context since it simultaneously explains the baryon asymmetry and the smallness of neutrino 
masses via the seesaw mechanism [3-8j. We add three heavy right-handed neutrinos Ni to the 
SM, which are assumed to have rather large Majorana masses Mi, close to the scale of some 
possibly underlying grand unified theory (GUT), Sgut ~ 10^^ GeV. In the early universe, the 
heavy neutrinos decay into leptons and Higgs bosons and create a lepton asymmetry, which is 
later on converted to a baryon asymmetry by the anomalous sphaleron processes [UdO]. The three 
Sakharov conditions [TT] that are necessary for a baryogenesis theory are fulfilled, that is lepton 
number L and B — L are violated, CP symmetry is violated in the decays and inverse decays and 
the interactions can be out of equilibrium. 
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Ever since the development of the theory 25 years ago, the calculations of leptogenesis dynamics 
have become more refined and many effects and scenarios that have initially been neglected have 
been considerecH. Notably the question how the hot and dense medium of SM particles influences 
leptogenesis dynamics has received increasing attention over the last years |13H20j . At high tem- 
perature, particles show a different behaviour than in vacuum due to their interaction with the 
medium: they acquire thermal masses, modified dispersion relations and modified helicity proper- 
ties. All these properties can be summed up by viewing the particles as thermal quasiparticles with 
different behaviour than their zero-temperature counterparts, much like the large zoo of single- 
particle and collective excitations that are known in high density situations in solid-state physics. 
At high temperature, notably fermions can occur in two distinct states with a positive or negative 
ratio of helicity over chirality and different dispersion relations than at zero temperature, where 
these dispersion relations do not break the chiral symmetry as a zero-temperature mass does. 

Thermal effects have been considered by references |13fl20j . Notably reference jl4j performs an 
extensive analysis of the effects of thermal masses that arise by resumming propagators using the 
hard thermal loop (HTL) resummation within thermal field theory (TFT). However, the authors 
approximated the two fermionic helicity modes with one simplified mode that behaves like a vacuum 
particle with its zero-temperature mass replaced by a thermal mas^. Due to their chiral nature, 
there are serious consequences to assigning a chirality breaking mass to fermions, hence the effects 
of abandoning this property should be examined. Moreover, it seems questionable to completely 
neglect the negative- helicity fermionic state which, according to TFT, will be populated at high 
temperature. We argue in this study that one should include the effect of the fermionic quasipar- 
ticles in leptogenesis calculations, since they behave differently from zero-temperature states with 
thermal masses, both conceptually and regarding their numerical infiuence on the CP-asymmetry. 
We calculate the full hard-thermal-loop (HTL) corrections to the CP-asymmetry in neutrino decays 
and, at higher temperature, Higgs boson decays, which have four different contributions, reflecting 
the four possibilities of combining the two helicity modes of the final-state lepton with the two modes 
of the lepton in the loop. As a comparison, we calculate the asymmetries for an approach where 
we approximate the lepton modes with ordinary zero-temperature states and modified masses, the 
thermal mass m^{T) and the asymptotic mass of the positive-helicity mode, ^/2mi{T). 

This paper is the first part of a two-paper series, where the second part is concerned with solving 
the Boltzmann equations using HTL-corrected rates and CP-asymmetries (22] . The present work 
deals with these corrections to the CP-asymmetries and is structured as follows: In section [21 we 
briefly review the imaginary time formalism of thermal field theory (TFT) and discuss the hard 
thermal loop (HTL) resummation. In section [3l we review our previous calculation for neutrino 
decays [23j and present a detailed analysis of the HTL-corrected rate for Higgs boson decays at 
high temperature. The CP-asymmetry for the different approaches is the main topic of section 21 
The CP-asymmetry in the two-mode approach consists of four different contributions due to the 
two possibilities for the leptons in the loops. We present some useful rules for performing calcula- 
tions with the fermionic modes and compare the analytical expressions for the CP- asymmetries in 
different cases. We restrict ourselves to the hierarchical limit where the mass of Ni is much smaller 
than the mass of N2, that is M2 ^ Mi and assume that the contribution of N^^ to the CP asym- 
metry is negligible. The temperature dependence of the CP-asymmetry is discussed in detail for 
the one-mode approach, the two-mode approach and the vacuum case. The differences between the 

^For an excellent review of the development in this field, we refer to reference [12| . 

^Moreover, an incorrect thermal factor for the CP-asymmetry was obtained, as has been pointed out in refer- 
ence [2T| . 
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asymmetries and the physical interpretation of certain features of the asymmetries are explained 
in detail. We summarise the main insights of this work in the conclusions and give an outlook 
on future work and prospects. In appendix \^ we derive frequency sums for the CP-asymmetry 
contributions. Analytical expressions for the CP-asymmetry are calculated in appendix [Bl while 
in appendix [O we present analytical expressions for the CP-asymmetry contributions of the two 
cuts through {N',i'} and {N' , which we did not consider in section HI since we are working in 
the hierarchical limit. We give an analytical approximation for the CP-asymmetry in Higgs boson 
decays at high temperature in the one-mode approach in appendix iDl 

2 Propagators at Finite Temperature 

When going to finite temperature [24j, one has to employ ensemble weighted expectation values of 
operators rather than the vacuum expectation values, so for an operator A we get 

{0\A\0)^{A)p^tT{pA). (2) 

There are two formalisms for calculating Green's functions at finite temperature, the imaginary 
time formalism and the real time formalism. Both are equivalent and we employ the imaginary 
time formalism, where the /co-integration is replaced by a sum over discrete energies, the so-called 
Matsubara frequencies. 

Naive perturbation theory at finite temperature can lead to serious conceptual problems, such 
as infrared divergent [25l[26] and gauge dependent [27tl28] results and results that are not complete 
to leading order. In order to cure these shortcomings, the hard thermal loop (HTL) resummation 
technique has been invented j291l30j. One distinguishes between hard momenta of order T and 
soft momenta of order gT, where g is the coupling constant of the corresponding theory. In a 
strict sense, this is only possible in the weak coupling limit where g <ti 1. If all external momenta 
are soft, then the bare thermal propagators have to be replaced by resummed propagators. The 
self-energies that are resummed are the HTL self-energies, for which all internal momenta are hard. 
For a scalar field with a HTL-self-energy 11, the resummed effective HTL-propagator A* follows 
from the Dyson-Schwinger equation in figure [J as 

iA* = iA + iA(-in)iA + ... 

= ' = ' (3) 

A-i - n K^-ml-W ^ ^ 

where A is the bare propagator, K the momentum and rriQ the zero-temperature mass of the scalar. 
The dispersion relation for this effective excitation is given by the pole of the propagator as 

k'^ = k'^+ml + U, (4) 

so we get an effective mass of rn^g = + nig where the thermal mass of the scalar is given 
by the self-energy, which is proportional to gT, rrig = 11 oc (gT)'^. It is possible to neglect the 
zero-temperature mass if ms ^ rriQ. 

For fermions with negligible or vanishing zero-temperature mass, the general expression for the 
self-energy in the rest frame of the thermal bath is given by [31] 

S(P) = -a(P)i*-6(P)^, (5) 

^We shamelessly stole our notation for the cuts in the vertex contribution from reference [2U| . 
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Figure 1: The resummed scalar propagator. 



where = (1, 0, 0, 0) is the four-velocity of the heat bath. The factors a and b are given by 

a{P) = ^ [tr(/'S) -potr(7oS)] , 
= ^ [^'tr (toS) - potr . 
In the HTL limit, the traces are given by [2 



(6) 



Ti = tr (i^S) =4:ml, 

T2 = tr (702^) = 2mp- In — 

P - P + 1 e 



(7) 



where the effective thermal fermion mass mp oc gT depends on the interaction that gives rise to 
the fermion self-energy. 

The resummed fermion propagator is then written as 



S*{K) 



^-ShtlW ■ 

It is convenient to rewrite this propagator in the helicity-eigenstate representation 
S*{K) = ^A+(i^)(7o - k • 7) + ^A_(i^)(7o + k • 7), 



(8) 



(9) 



where k = k/A;, and 



A±(K) 



-ko±k + 



m. 



k 



±1 



ibfco — k kQ + k 
■ In 



2k 



kQ — k 



(10) 



This propagator has two poles, the zeros of the two denominators A±. The poles can be seen 
as the dispersion relations of single-particle excitations of the fermions that interact with the hot 
plasma. 



ko = 0±ik) . 



(11) 



We have presented an analytical expression for the two dispersion relations making use of the 
Lambert W function in reference |23] . The dispersion relations are shown in figure El 

Note that even though the dispersion relations resemble the behaviour of massive particles and 
= mp for zero momentum k, the propagator S*{K) ([9]) does not break chiral invariance like a 
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Figure 2: The two dispersion laws for fermionic excitations compared to the standard dispersion 
relation 0^ = /c^ + rnj^. 



conventional mass term. Both the self energy '^{K) ([5]) and the propagator S*{K) anticommute 
with 75. The Dirac spinors that are associated with the pole at ko = 0+ are eigenstates of the 
operator (70 — k • 7) and they have a positive ratio of helicity over chirality, = + 1. The spinors 
associated with ko = 0_, on the other hand, are eigenstates of (70 + k • 7) and have a negative 
helicity-over-chirality ratio, = — 1. At zero temperature, fermions have = + 1. The introduction 
of a thermal bath gives rise to fermionic modes which have = — 1. These modes have been called 
plasminos since they are new fermionic excitations of the plasma and have first been noted in 
references [31113^ . 

We can introduce a spectral representation for the two parts of the fermion propagator [35] , 



00 



A.(K)=/ d0-J^^^, (12) 



00 



where the spectral density p±{0,k) |32y36j has two contributions, one from the poles, 

p^°''(0, k) = Z±(0, k) (0 - 0±{k)) + Z^{0, k) (0 + 0^{k)) , (13) 
and one discontinuous part. 



9f'{0, k) = \ml{k^0)^ _ ^ ^^^2 _ 

{k{0^k)-m?p[QQ{x)^Qi{x)]] + [i7rm2,(l^x)] 



where x = 0/k, 9{x) is the heaviside function and and Qi are Legendre functions of the second 
kind, 

1 X + 1 

<5o(2;) = -ln -, Qi{x) = x Qq{x) - I . (15) 

2 X — \ 



The residues of the quasi-particle poles are given by 

0i{k)- e 

2 m 



Z±{0, k) = ^\ ' ^ , where Z+ + Z_ = 1 . (16) 
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Figure 3: The momentum-dependent effective masses m±. 

One can describe the non-standard dispersion relations 0± by momentum-dependent effective 
masses m±{k) which are given by 

m±{k) = ^0l{k)-k^ = yj2Z{0,k)mF . (17) 

These masses are shown in figure El 

Considering gauge theories, one might also have to use HTL-corrected effective vertices that 
are related to the propagators by Ward identities [M|- We do not consider these vertices since we 
are only looking at Yukawa vertices. In the HTL framework, it is sufficient to use bare propagators 
if at least one of the external legs is hard. However, it is always possible to resum self-energies 
and thus capture effects which arise from higher-order loop diagrams and take into account the 
appearance of thermal masses and modified dispersion relations in a medium. In fact, since the 
effective masses we encounter do typically not satisfy the condition mgff <C T but are rather in the 
range nies/T ~ 0.1 - 1, the effect of resummed propagators is noticeable even when some or all 
external momenta are hard. In summary, we always resum the propagators of particles that are 
in equilibrium with the thermal bath, which are the Higgs bosons and the leptons in our case, in 
order to capture the effects of thermal masses, modified dispersion relation and modified helicity 
structures. This approach is justified a posteriori by the sizeable corrections it reveals, similar to 
the treatment of meson correlation fuctions in reference |37j . 

In leptogenesis, the leptons and Higgs bosons acquire thermal masses that have been calculated 
in references [HHIMIEEIEH] and are given by 

The couplings denote the SU(2) coupling (72, the U(l) coupling gy, the top Yukawa coupling yt and 
the Higgs self-coupling A, where we assume a Higgs mass of about 115 GeV. The other Yukawa 
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Figure 4: N decay via the optical theorem with dressed propagators denoted by a blob. 



couplings can be neglected since they are much smaller than unity and the remaining couplings are 
renormalised at the first Matsubara mode, 27rT, as explained in reference [\Al and in reference |4Uj 
in more detail. The zero-temperature Higgs boson mass is negligible compared to the thermal mass 
and the SM fermions do not acquire a zero-temperature mass since the temperature is above the 
electroweak symmetry breaking scale. The heavy neutrinos A^i do acquire a thermal mass, but since 
the Yukawa couplings are much smaller than unity, this effective mass can be neglected compared 
to the zero-temperature mass. 



3 Neutrino and Higgs Boson Decays 

In leptogenesis, we add three heavy right-handed neutrinos Ni to the SM, which are assumed to 
have large Majorana masses Mi close to the GUT scale, -Bqut ~ 10^^ GeV. The additional terms 
in the Lagrangian are 

C = i m^jf^Ni - liaN^ireabO -^Yl MiNi^f + h.c. , (19) 

i 

where the Higgs doublet cp is normalised such that its vacuum expectation value (vev) in 

= (°) (20) 

is ~ 174 GeV and ka is the Yukawa coupling connecting the Higgs doublet, the lepton doublet 
and the heavy neutrino singlet. The indices a and b denote doublet indices and eab is the two- 
dimensional total antisymmetric tensor that ensures antisymmetric S'C/(2)-contraction. 

We have discussed the HTL corrections to neutrino decays A^i — )• HL in detail in reference [23] . 
When the temperature is so high that > Mi, the neutrino decay is kinematically forbidden 
in the HTL-approximatior|^, but the decay of Higgs bosons into neutrinos and leptons becomes 
possibl^. The rate for the Higgs boson decays can be calculated in the same way as the rate 
for the neutrino decays, by cutting the A^i-self-energy with resummed lepton and Higgs boson 
propagators in figure [H According to finite-temperature cutting rules j4HI42]. the interaction rate 
reads 

r(P) = --^ tT[{^ + Ml) Im S(po + ie, p)]- (21) 

^It has been shown in reference [16], that the decay is still allowed if one considers the effect of coUinear external 
momenta. 

'^The lepton decay is not possible, since > mi for all temperatures 
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At finite temperature, the self-energy for a neutrino with momentum P is given by the Matsubara 
sum 



E(P) = -4(/tOnT 

ko=i {2n+l)nT 



(2^ 



PLS*{K)PnD*{Q), 



(22) 



where S* and D* are the HTL-resummed lepton and Higgs boson propagators in equations ([3]) 
and ([9]), Pl and Pr are the projection operators on left- and right-handed states, K is the lepton 
momentum and Q = P—K the Higgs boson momentum. We have summed over the two components 
of the lepton and Higgs doublets, over particles and antiparticles and the three lepton flavours, so 
we are looking at the processes H o NiL^ where the notation H and L indicates that we are 
considering both I and I. 

Since the leptonic quasi-particles are the final states, we are only interested in the pole contri- 
bution of the lepton propagator and we get for the Matsubara sum [23J 



feo 



1 
2a;„ 



+ 



2m] 

uj^ikf - e 

2mj 



l + /^(0g) -/K0±) , U{0q) + fe{0±) 



PO - UJ±{k) - UJq po - U}±{k) + UJq _ 

Ui0q) + fei0T) ^ 1 + f4'i0q) - M0t) 



PO+UJ^{k)-UJq po + UJ^{k) + UJq 



(23) 



where ujq 



+ m^, is the energy of the Higgs boson, 0±{k) denotes the two lepton dispersion 



relations, f(),{0q) = [exp(0q) — 1]~ is the Bose-Einstein-distribution for the Higgs bosons, fe{0±) = 
[exp(0-|-) -|- the Fermi-Dirac distribution for the leptons and =1/T. 

The four terms in equation (I23p correspond to the processes with the energy relations indicated 
in the denominator, i.e. the decay A'^i — )• HL, the production NiH — )• L, the production NiL — )• H 
and the production of NiLH from the vacuum, as well as the four inverse reactions [41]. We are 
only interested in the process H -f-)- NiL, where the decay and inverse decay are illustrated by the 
statistical factors 



U + fi = Ui^-h) + ii + U)fe. 



(24) 



The decay is weighted by the factor /^(l — fi) for absorption of a Higgs boson from the thermal 
bath and induced emission a lepton, while the inverse decay is weighted by the factor (1 -|- f^)f£ 
for induced emission of a Higgs boson and absorption of a lepton from the thermal bath. Our term 
reads 



ko 



H^NiL 



1 

2uJn 



2mi 



PO + LO-h 



(25) 



where /i = ±1 denotes the helicity-over-chirality ratio of the final-state leptons. The angle r] 
between the final-state neutrino and lepton is given bjH 



Vh 



Wp + 



(26) 



*^Note that the physical three-momenta of the initial-state Higgs boson and the final-state neutrino are — q and 
-p since we were starting from the neutrino self-energy and not from the Higgs boson self-energy. 
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where 



= m2 - m2 - {ul - k') . (27) 

In order to clarify the momentum relations, we revert the direction of the three- momenta q — )• — q 
and p — )• — p so that they correspond to the physical momenta of the incoming Higgs boson and 
outgoing neutrino. The matrix element can then be derived as 

\Mh{P, K)f = 4. {lU)nZhOJh{Po -hp-k) = A {lU)nP^K > (28) 

where we have introduced a chirally invariant four-momentum Kj^ = Z/j0/i(l, /ik) for the lepton. 
This matrix element looks the same as the matrix element for neutrino decays in reference |23j . 
since the momentum flip p — )• — p compensates the helicity flip h ^ —h and 

is the residue of the modes and the angle for the reverted physical momenta reads 

Vh = ^ i'^Pou^h - . (30) 

We have derived the matrix element for the Higgs boson decays starting from the neutrino 
self-energy since the neutrino spinors are not affected by thermal corrections, so we could extract 
the matrix element from the expression for the neutrino interaction rate in equation (|2ip . It 
is also possible to derive the matrix element from the Higgs boson self energy, even though the 
Higgs bosons are affected by thermal corrections, but their external states are the same as the 
vacuum states since the thermal propagator in equation ([3]) has the same structure as the vacuum 
propagator with a different effective mass. 

For leptons, the situation is different, because the structure of the HTL-propagator in equa- 
tion Q is structurally different from the vacuum propagator. This means that the external spinors 
will have a different structure due to the modified helicity properties of the HTL propagator. We 
derive important properties of the effective lepton spinors in section [4Tl 

Integrating over all neutrino momenta p, the decay densities for neutrino and Higgs boson decay 
are given by 

jiN^HL) = J dpNdpHdpU27r)^S\pN-PL-PH)\MH\'f^\l-fr){l + rj?) (31) 



-f{H^NL) = J dpNdpHdpLi2n)^6\pN+PL-PH)\MH\H'^-fT^)a-fr)fi?, (32) 



and 



where dp = d^p/[{27r)^2E] and the matrix element is defined given by equation (128p . 

In figure [5l we compare our consistent HTL calculation to the one- mode approximation adopted 
by reference [H], while we add quantum-statistical distribution functions to their calculation, which 
equals the approach of using an approximated lepton propagator — mi) [l3]. We have shown 

the decay density for the Higgs boson decays in reference 01], but present a much more detailed 
analysis in this work. In addition, we show the one-mode approach for the asymptotic mass \f2rag. 
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T/jTljv 

Figure 5: The decay densities for the neutrino and the Higgs boson decay. We show the one-mode 
approach with the thermal mass as 7^,^ and with the asymptotic mass as 7^^^; Also the T = 
rate 70 and our two modes 7-1- . The temperature thresholds are explained in the text. 

We evaluate the decay rates for Mi = 10^" GeV and normalise the rates by the effective neutrino 
mass mi = {lH)iiv'^/Mi, where v = 17 A GeV is the vacuum expectation value of the Higgs field. 
This effective mass is often taken as rhi = 0.06 eV, inspired by the mass scale of the atmospheric 
mass splitting. 

In the one-mode approach, the decay is forbidden when the thermal masses of Higgs boson and 
lepton become larger than the neutrino mass, Mi < + or Mi < \plm( -|- m^. Considering 
two modes, the kinematics exhibit a more interesting behavior. For the plus-mode, the phase space 
is reduced due to the larger quasi-mass, and at Mi = \/^m£ + m^, the decay is only possible 
into leptons with small momenta, thus the rate drops dramatically. The decay into the negative, 
quasi-massless mode is suppressed since its residue is much smaller than the one of the plus-mode. 
However, the decay is possible up to Mi = m^. Due to the various effects, the two-mode rate differs 
from the one-mode approach by more than one order of magnitude in the interesting temperature 
regime oi z = T /Mi > 1. The \/2 m^-calculation is a better approximation to the plus-mode, but 
still overestimates the rate, which is due to the different structure of the matrix elements in the 
one-mode and the two-mode approach, that is the helicity structure of the quasiparticles. The 
residue also reduces the plus-rate, but the effect is smaller since is usually close to one. 

At higher temperatures, when > Mi + m±{k), the Higgs can decay into neutrino and lepton 
modes and this process acts as a production mechanism for neutrinos [14]. The decay (p — )• Nl^ 
is possible when > Mi, while the decay into is possible when > Mi + m£. As for 
low temperature, the rate 74. is unsuppressed only when > Mi + \f2m,i. Our decay density 
approaches the decay density of reference [H] at high temperatures, but is about a factor two below. 
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This can be explained by the fact that the phase space is smaller due to the larger mass of the 
lepton, m( < mh{k) < \f2m(,. Again, the asymptotic mass calculation is a better approximation 
but still gives a larger rate due to the matrix element and, to less extent, the residue. We see that 
the decay rate rises as T'^, instead of ^ as for the vacuum rate 70. In the vacuum calculation, 
the squared matrix element is proportional to M\. In the finite temperature calculation, it is 
proportional to = — mf^{k) — Ml, so the dominant contribution is proportional to m| ~ 
and the rate rises by a factor faster than the vacuum rate 70. 

Summarising, we can distinguish five different thresholds for the thermal decay rates we dis- 
cussed. Going from low temperature to high temperature, these are given by the following condi- 
tions: 





Ml 


= V2mi + 


rpN 


Ml 


= me + m^, 




Ml 


= , 


rp<t> 

-'o 




= + Ml , 






= V2me + Mi 



4 CP- Asymmetries 

4.1 Defining CP-asymmetries at finite temperature 

Let us turn to calculating the C!P-asymmetry in Ni decays. We denote the decaying Ni by N 
and the N2 in the loop by N' and we assume that the contribution of to the CP asymmetry is 
negligible. At T = 0, the CP-asymmetry is defined as 

_ r(jV ^ cPi) - T{N ^ 

where T are the decay rates of the heavy A'^s into Higgs boson and lepton doublet and their 
CP-conjugated processes. At finite temperature, we have to calculate the CP-asymmetry via the 
integrated decay rates, 

^ 7^>°(A^04)-7^>°(A^^'A4) 
^'^^ > ^T>o(^r^04) + ^T>o(^r^<^4)' ^^^^ 

where we define the CP-asymmetry for each lepton mode, denoted by h. We have 

7^>° = / ^,fN{PN)r^>'{P^), (36) 

where /jv is the distribution function of the neutrinos and the neutrino momentum. In the 
zero-temperature approximation, we write 

Ml 

r(pM) = ^Frf, (37) 
Po 

where Mi and pQ are the mass and the energy of the neutrino and Fj-f is the decay rate in the rest 
frame of the neutrino. The integration over the momentum cancels out and the CP asymmetry 
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via 7 is the same as via T. At finite temperature, however, the thermal bath breaks Lorentz 
invariance and the preferred frame of reference for calculations is the rest frame of the thermal 
bath. The momentum dependence of the decay rate cannot be formulated as in equation ()37p and 
the CP-asymmetry as defined in equation (j34p is momentum dependent, therefore the definition in 
equation ([35|) is the appropriate one. 

The CP asymmetry in equilibrium can be written as 



J J^ppiJNK^ Dh + i-Dh 



(38) 



where T^h = V{N — l^'P) ^'^^ ^Dh = r(A^ — )■ Ihcj)) are the decay rate and the CP-conjugated 
decay rate. 

The decay density is written as 



7Dh 



where 



Zd = {1- /7v)(1 + a - fi) = (1 + /0)(1 - fi) 



(39) 



(40) 



is the statistical factor for the decay, with Bose-enhancement and Fermi-blocking. In the denomina- 

2 2 

tor of the CP-asymmetry, it is sufficient to take the tree- level matrix element. 
The CP- asymmetry reads 





2 




A^trce 




A^tree 



fdE dk^ In Zn{\Mh\^-\Mh\'') 
2 JdEdk^ fN ZD\Mh\^ 

7.(iV^L^) 4(2^)3/^^^^^^^ 



Ml 



Ml 



(41) 



The CP-asymmetry arises as the interference between tree-level and one-loop diagrams in the 
decay, so we write M = A^o + Mi^ where A^o is the tree-level amplitude and Mi the sum of 
all one-loop amplitudes. The matrix elements can be decomposed as Mi = Xili such that the 
CP-conjugated matrix element is Mi = Here, Aj includes the couplings and Ij accounts for 

the kinematics. Thus, 



where Xcp = AqA* and Icp 



lA^P 



lA^p 



-4 Im Xcp Imlcp, 



(42) 



4.2 The vertex contribution 

Calculating the imaginary part of the kinematic term Imlcp amounts to calculating the imaginary 
part of the one-loop diagram since the tree-level diagram is real. There are two one-loop diagrams 
for the neutrino decay, the vertex diagram and the self-energy diagram. The vertex diagram is 
shown in figure El along with the momentum assignments. The coupling is 

Xcp = XoX\ = [{X'^X),k?9su{2). (43) 
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Figure 6: The momentum assignments for the vertex contribution to the CP asymmetry. The sohd 
Hnes without arrows are neutrinos, the ones with arrows the leptons and the dashed hues the Higgs 
bosons. Ah momenta are flowing from left to right and P' as indicated. 



where gsu{2) = 2 denotes the sum over the Higgs and lepton doublets, j = 1 is the decaying neutrino 
family, A; = 2 is the family of the neutrino in the loop and we have summed over all fermion spins 
and the lepton families, both external and in the loop. Moreover, 

= -iy [M2AN,A^,{uePRUN){uNPRSePLUi)]* , (44) 

where p', q' and k' are the neutrino, Higgs boson and lepton momentum in the loop, A^/ = 
(P'^ — MI)"^ is the denominator of the loop neutrino propagator, A^i accordingly for the loop 
Higgs, Si' is the loop lepton propagator and un and U£ are the external neutrino and lepton spinors. 

The external fermions are thermal quasiparticles and can be written as spinors [24j which are 
eigenstates of (70 =F k • 7) and have modified dispersion relations. We have shown in reference |23] 
that 

i ^ \MUP,Kr = {PoTPV±) , (45) 

s t 

where s denotes the spin of the neutrino. We can also write the matrix element as 

\Y.\^±iP^Kt = lY.9'i4PR^N){u%PLuf). (46) 

s s 

From equations (j45p and ()46p we derive a rule for multiplying the spinors of the lepton states, 

uf {K)uf (K) = Z±oj± (70 ^ k • 7) , (47) 

where 

Z± = (48) 

is the quasiparticle residuum. For the antiparticle spinors v, we replace K by —K and get 

vf{K)vf{K) = -Z±a;±(7o ± k • 7) • (49) 
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The HTL lepton propagator is given in equations ([9]) and (|T0|) and the Higgs boson propagator 
in equation ([3]). At finite temperature, we sum over the Matsubara modes, 



"-0 

where 



ft,,. 



fco = (2n + l)7riT, (51) 

since we are integrating over a fermion momentum. 
The spin and helicity sum are evaluated as 

Y.iu'lPRU^W^PRS^: Pl4) = -Y1 ZhOOhMiAh'il - hh'k ■ k'), (52) 

s,h' h' 

where h and h' are the ratios of hehcity over chirahty for the external and the loop lepton. The 
integral reads 

M2M^Zh0h [A^, A^, a;,] * , (53) 



k'Q,h' ' 



where H± = l± hh'kk'. 

The frequency sum is calculated in detail in appendix IA.2I We are only interested in the 
contribution from the pole part of the lepton propagator and the explicit expression is 

k'g h' 



+ 



[< - <) 4 - - <) A', - (< - B%) ^? + {Bl - B%) Af] H^} , 

(54) 

where the factors -Bat/^ and Ai are given by 

=Pq-Uj' -UJqt, Nf = ko- U)q, - UJpt, = Qq - Uj' - UJp, , 

N% = Pq + J + UJqi, = ko + Uq' + Up/ , = Qq + Uj' + UJp' , 

= PO - ^' + ^q', Nf =ko- UJq> + Up,, = - uj' + UJp/ , 

Nff = Po + Uj' - UJq/ , N^' = ko+UJq' - UJp>, N^' = + J - UJp, , (56) 

= 1 - Z^' + /<^' ) = 1 ~ ff-' ~ fN' , 

= -(1 - fi' + f<t>') ' = ~(1 ~ f(-' ~ fN') , 

= ~ife' + /</-') ) = -{fe' - In') , 

zff = fe' + f<j,' , = fii - f^i . (57) 
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Figure 7: The momentum assignments for the self-energy contribution. The sohd Unes without 
arrows are neutrinos, the ones with arrows the leptons and the dashed lines the Higgs bosons. 



4.3 The self-energy contribution 

For the self-energy contribution, the integral Ig = ly is the same as for the vertex contribution, only 
the momentum relations are different (cf . figure [7]) . The left diagram does not give a contribution 
since the combination of couplings, |(A^A)jfcp, does not have an imaginary part. 

In order to carry out the Matsubara sum, we use the Saclay-representation for the propagators. 
For the Higgs propagator it is given by 

A'^ = - f dr e-^o^-^lll + /<^,K,)]e-"''" + /^^KOe""'"}, (58) 



where ojqi = y ^'^ + '^n^ is the on-shell Higgs energy with the thermal Higgs mass and f(f,i is 
the Bose-Einstein distribution for the Higgs bosons with energy ojqi. For the lepton propagator the 
Saclay representation is given by equations (pl|) and (j95|) . The neutrino propagator simply reads 

since the internal neutrino momentum P' is the same as the external neutrino momentum P. As 
usual, we can write = i (2m -|- l)7rT as Matsubara frequency and later on continue it analytically 
to real values of po- Iii particular e^"'^ = — 1. 
We can calculate the frequency sum directly. 



T e^o-e^o"' = e^ »"5(r' - r) (60) 

k' 

and get 

/OO 1 
^du:'p'(h')^[B^-B';) . (61) 

Alternatively, we can use equation (199p and write 

^ 2^ K',s (ko, ^ ) As^, {po - ko, ujg,) = Zm— ; ^ . (62) 

~f 2cjq/ Pq - SUJsh' - S^iOJqi 

Kg 

Both calculations lead to 

A^,A'//_ = ^Z,,[{B^ - Bi)H^ + {Bf; - B%)H^] . (63) 

k'g h' h' ^'i' 
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Figure 8: The cuts through the vertex contribution at finite temperature. The cuts are closed to 
form circles and the line that denotes the decaying particle in the corresponding 1 — )• 2 process is 
indicated by a blob. 



4.4 Imaginary parts 

The terms B'^^^ and Af in the vertex contribution in equation (I114p correspond to the three 
vertices where the denominator fulfills certain momentum relations when set to zero: the Bjsf- 
terms correspond to the vertex with an incoming A'^i and {£' , (j)'} in the loop, the S^-terms to the 
vertex with an outgoing (p and {N2,i'} in the loop, and the ^^-terms to the vertex with an outgoing 
i and {A''2, (j)'} in the loop. As an example, the term 

= ^-f^' + f^' (64) 

PO -UJ' - idq' 

corresponds to the incoming neutrino decaying into the lepton and Higgs boson in the loop. Thus, 
the terms correspond to cuttings through the two loop lines adjacent to the vertex, however, a 
correspondence with the circlings of the RTF [21^i24j is not obvious. Among these cuts, only the 
ones which correspond to a A^i or decaying into a Higgs boson and a lepton are kinematically 
possible at the temperatures where neutrino decay is allowed, that is where Mi < m^. These terms 
are B§, Af and Bf. 



Regarding the cuts 

The diagrams develop an imaginary part when one of the denominators of the relevant terms B^, 
and B^ vanishes. The contributions from these denominators, , and A'^^ , correspond 
to the three possible cuts shown in figure [H The contribution from is the only possible cut 
at zero temperature. At finite temperature, the other two cuts correspond to exchanging energy 
with the heat bath. When choosing the imaginary parts corresponding to these two cuts, the loop 
momentum K' is of the order of M2. Since we assume a strong hierarchy M2 ^ Mi, the thermal 
factors fcf,',fi' and /^r/ are suppressed by the large loop momentum and the contributions become 
very small. In fact, they turn out to be numerically irrelevant in the hierarchical limit. The physical 
interpretation of this is as follows: Consider for example the cut through {^', which is given 
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by a vanishing denominator A'^^'. The corresponding thermal weighting factor is the numerator 

— fi' ~ = /^'(l ~ fN2) ~ (1 ~ fl')fN2- It corresponds to two processes: absorption of a 
neutrino from the thermal bath and induced emission of a lepton, or absorption of a lepton and 
induced emission of a neutrino. The phase space distribution of the A'^2S in the bath is suppressed 
due to their large mass and also the distribution of ^s that have momenta large enough to fulfill 
momentum conservation in the process is suppressed, so the process is suppressed. Therefore, the 
thermal factors suppress the contribution from the A'^2-cuts. Only when we have degenerate masses 
M2 > Ml, these cuts will give a contribution similar to the one from . In this cas^, the energy 
and temperature scales that correspond to A''! and N2 processes are not clearly separated and one 
has to account for the possibility of an asymmetry creation by A''2 as well. Implications of these 
cuts were discussed in reference [20]. We do not consider the influence of this cuts, since we are 
working in the hierarchical limit, but we present the analytical expression in appendix [Cl 



Vertex cut through (f)'} 

The imaginary part from , which implies cutting through the lepton and Higgs boson in the loop, 
is the only cut that is also possible at zero temperature and the only vertex cut that contributes 
in the hierarchical limi10. We denote the angle between p and k' with 77', 



Then 



Im (y"^ dr?'^) = -vr dV<5« ) = -vr ^6{v' - v'o) (66) 



where the angle is 



and 



pk' ' 



We get 



^ ^fdk'dcp'^Zh'Z^{Ai-Af)H.. (69) 



167r2 ^ 
h' ■ 



''Note that there is a mass range for M2 where we have a contribution from the cuts but no resonant enhance- 
ment by the self-energy contribution, which becomes relevant when AM = M2 — Mi <C Mi. This mass range is at 
AM ~ Ml > r 

^'^The corresponding CP-asymmetry has been calculated in reference [14], but with a thermal factor 1 — fii + f^r — 
2fiif^i instead of the correct 1 — + f^i. For details, see reference |21| . 



17 



It is sufficient to perform the integration over (j)' from to tt since cos (j)' in is the only quantity 
that depends on cf)' . 

We note that we can write 

4 - Af = ^ = 2^,,A]^^ , (70) 

(/CO-Wg') -Wp' 

where can be viewed as the propagator of the internal neutrino, since we can interprete 

the contribution we are looking at as putting the internal Higgs boson on-shell and thus we have 
ko - ujq' =ko-q'Q= p'q. 

The analytic expression for the CP-asymmetry as defined in equation (j35|) is worked out in 
app endix IB . 1 1 and given by 

Im{[(AtA)dnMiM,E.Jdi^dMfc7;d.^-fcF^Z,^Z^Z,,(^r<)g- 
"''^ ^ 5c(AtA)ii 47r2 jAEdkkfNZDZh{po-hpr^) ' ^ 

where gc = indicates that we sum over N ^ (pi and N ^ and F/v/i = /^(l + /^^)(1 — fg^) 
the statistical factor for the decay. 

Self-energy cut 

For the self-energy diagram, only contributes. Taking 7]' as the angle between p and k', we get 



Comparing this expression with the contribution from A^^v in equation ()69p . we see that calculating 
the self-energy contribution amounts to replacing A]^i^ by = (M^ — M"^)^^ in the A^Tv-vertex 
contribution. If M2 ^ Mi, we get 



SO the self-energy contribution is twice as large as the vertex contribution, eg ~ 2ev, where the 
factor two comes from the fact that we have two possibilities for the components of the SU{2) 
doublets in the loop of the self-energy diagram. This resembles the situation in vacuum. 

The analytic expression for the CP-asymmetry is worked out in app endix IB . 2 1 and given by 

Im{[(AtA)i2]n M^M2 1 ZH'I'iEdkdk'£d<t>'kF^^lZ^!^Z^Z,,H. 
^""^ ' <?c(AtA)ii M2-M|27r2 j dEdkkfNZoZnipo - hpr^) ' ^ ^ 

Symmetry under lepton-mode exchange 

We can use equation (j40p and collect all factors that depend on k and k', 

(1 + /</, - + /</,' - fi')ZhZt,,kk'Al'^H_ , (75) 
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Figure 9: The product of diagrams that corresponds to the vertex contribution of the CP asymmetry 
at low temperature. It is symmetric under the exchange of the leptons I and (! . 

where we have suppressed the indices for hehcity-over-chirahty ratios h and h! . The internal 
neutrino momentum p = k + k' — p' is symmetric under a replacement of k and k' and likewise 
the difference w — lOqi = a; + w' — po- The Higgs boson momenta q = p — k and q' = p — k' are 
also exchanged when we exchange k and k'. Thus, the C!P- asymmetry for the vertex contribution 
is symmetric under an exchange of the internal and the external lepton. This can be understood 
as follows: Taking the imaginary part of ^AQ^A\ by putting the internal lepton and Higgs boson 
on-shell corresponds to calculating the product of the amplitudes of two decays and one AL = 2 
scattering with a neutrino in the u-channel, as shown in figure [9l It can easily be checked that this 
symmetry also holds for the self-energy diagram, where the corresponding AL = 2 scattering has 
a neutrino in the s-channel. 



4.5 The CP- asymmetry at high temperature 



At high temperature, where we have the decays of Higgs bosons, the CP-asymmetry on amplitude 
level is defined as 



(76) 







\M{^ ^ iV4) 





The external momenta are now related as go = po + ^o- The momentum assignments are shown in 
figure [ini We take q and p as the three-momenta of the initial-state Higgs boson and the final- 
state neutrino as in section [3l this way we can directly use the results from the CP-asymmetry in 
neutrino decays. The matrix elements are the same as for the low temperature case, so A^(</> — ?■ ISSi^) 
corresponds to M.{N — )• ^Ih)^ just the energy relations are different. The self-energy contribution 
from the external neutrino line is the only CP-asymmetric self-energy, the other self energies do 
not exhibit an imaginary part in the combination of the couplings. The couplings read 



Im 



{A^f} 



9SU{2) 



Im 



AtA 



21 



-9SU{2) 



Im 



AtA 



12 



(77) 
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Figure 10: The vertex and the self-energy contribution for the (j) decay. 



The integrals for the vertex and the self-energy contribution are 



T<P t9' 

-'o-'i 



^ J (2vr)' 



(78) 



where we remember that A^r/ = l/{Mf — M|) for the self-energy graph. 
The frequency sum for the vertex diagram reads 

fc(', h' 



+ 



}■ 

(79) 



Since we have M2 ^ Mi, we also have M2 ^ m,^ in the relevant temperature range, so the 
possible contributions are from A''^ , N^' and A^^l^. Again, the N2 cuts can be neglected because 
they are kinematically suppressed. When taking the discontinuity of the diagrams, we get for the 
angle between p and k'. 



where 

so we arrive at 

/TV N\ ImXcp 



(80) 



(81) 



■yhlehjy 



4(27r)^ 



M1M2 J2 I dEdkdk' I " d(l)'kF^hZh-—Zh'Zf^{Af - A'DH^ , (82 



^^If niA, 3> A/2, we would have contributions from Nf and Nt instead 
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where we can write 

< - = ^2 = 2-p'A];f . (83) 

{Kq +U}qi) — 

Contrary to the CP-asymmetry in neutrino decays, this expression can not strictly be seen as the 
propagator of the neutrino since the contribution does not correspond to a zero temperature cut but 
is a pure thermal effect induced by the presence of leptons and Higgs bosons in the thermal bath. 
This is illustrated by the factor = f^i + fi' = /<^'(1 — /^') + (1 + AO/^'El; which describes the 
absorption of a Higgs boson and the stimulated emission of a lepton and the opposite process, the 
absorption of a lepton and the stimulated emission of a Higgs boson. Compared to low temperature, 
we have replaced A^f'Z^ by A^fzff. 

For the self-energy diagram, the frequency sum is given by 

E A7v'A^.A'i/_ =Ar,~Yl ^y^' (^^ - ^A^') + ^+ (^^ -^n)\^ (84) 

fc^ h' h' 

after taking the discontinuity, the CP-asymmetry reads 



= / dEdMk' J\4>'kF,,Z,h'Z,,{-Zt;)H_ , (85) 



where 

F,, = f;''{i-Q){i-r^). (86) 

Compared to low temperature, we have replaced Z^ by zfj . The self-energy contribution is given 
by replacing A^f by A^f = [Mf — M|) in the vertex case. For M2 3> Mi, we have 

so the relation ~ 2€v also holds for the Higgs boson decays. 

Using /^(l — fi) = [f^ + fi)fN, the terms that depend on the lepton momenta k and k' are 



iU + fi)iU + h')ZhZh'kk'A%'"'H_ . 

where now p' = k + k' -|- p and uj + ooqi = lo + uj' + so the CP-asymmetry in Higgs boson decays 
is symmetric under exchanging the internal and external lepton as well. 

4.6 One-mode approach 

We also calculate the CP-asymmetry within the one-mode approach where we treat the thermal 
mass like a kinematical mass and use lepton propagators {^—m()^^ or {^ — ^/2m()~^ as in section[3l 
The spin sum corresponding to equation ()52p then reads 



Y,{u}PRU%){u%PRSePLu}) = IMiAeRf^K'^ , (89) 



^^Reference |14) obtains a different factor /^z — fii ~ 2/^/ /^/ due to an incorrect choice of cutting rules as explained 
in reference [21] , 
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where A/ 



2 \-l 



pole 



. In the frequency sums in equations (|102p and ([62]) . we replace A^^, 



by the usual decomposition A<j^£/, which means replacing Zgh' by —s/{2lo') on the right-hand sides. 
One can check that in the final expression for the CP-asymmetry, this amounts to replacing the 
sum of the helicity contributions 



hh' 



hh'O by ^ 

LOkiOk' 



1 



kk' 



(90) 



This means that in the two mode treatment, it is forbidden for the external and internal lepton to 
be scattered strictly in the same direction if they have the same helicity or in the opposite direction 
if they have opposite helicity. For the one-mode approximation this is not the case since oJkOJk' is 
always larger than kk' . This result illustrates that the leptonic quasiparticles still behave as if they 
are massless in terms of the helicity structure of their interactions, while the one-mode approach 
is not able to describe this behaviour. 

For the CP-asymmetries in the decay densities we get 



(At 



m)v 



ImAcp 



7" 



(27r)^ 



{N-. 

d-Ed/cd/c' 

d£;dA;d/c' 



-— — r^MiMaV 
2(27r)5 

ImAcp , ^ , ^ \ - 
^MiM2 

h' 



TP 



K-K' 



LJk^k' 



^'^J^kk'Alf^^ 

P COk^k' 



ImXcp 
'72^ 



f M1M2 ^ / dEdkdk' / d^' 



p 



K-K' 

UJk^k' 



(91) 



where F^h = fN{^ - fN){l + /</.- fih), F^h = /Af(l - fN)iU + hh)- 

We can examine the high temperature behaviour of the one-mode approach by calculating the 
CP-asymmetry in the matrix elements of a Higgs boson at rest, where we assume that Mi, mi ^ 
m0 ^ M2. The algebra is worked out in appendix |Dj The result reads 



Q J. 



leo, 



(92) 



where eo is the CP-asymmetry in vacuum and = m^/T. Assuming that (7^ <C 1, we get 



32 

72 



(93) 



Taking g^p = m^/T ~ 0.42 for T = 10^^ GeV and using the more accurate term in equation (j92p . 
we get e/eo ~ 70, while we get e/eo ~ 90 for equation We view this result as a rough 

approximation of the value of the CP-asymmetry in Higgs boson decays at high temperature. Both 
our approximation and the numerical solution of the exact expression in the next section give a 
factor of 100 difference to the CP-asymmetry in vacuum. 
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4.7 Temperature dependence of the CP-asymmetries 



We show the temperature dependence of the CP-asymmetries in neutrino decays in the full HTL 
calculation and in the one-mode approach for m£ and \f2mi in figure[TTJ We choose M\ = 10^*^ GeV 
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Figure 11: The CP-asymmetries in neutrino decays normalised by the CP-asymmetry in vacuum 
and the total decay density in vacuum, A7/(7Q°*eo). We choose Mi = 10-'^° GeV and M2 » Mi. 
The term A'yh-^h^ denotes the difference between the decay rate and its CP conjugated rate, which 
is proportional to the CP-asymmetry. Here, hi denotes the mode of the external lepton, while /12 

denotes the mode of the lepton in the loop. For example, A7_| = 7(A'' — )• — 7(A'' — )• (^^+), 

where a minus-mode lepton is present in the loop, /^^rrn and ^7^^^ denote the rate differences 
for the one- mode approach with a thermal mass and an asymptotic thermal mass \f2mi. 

and normalise the asymmetries by the product of the CP-asymmetry at zero temperature and the 
total decay density in vacuum, eo7o°* ■ discussed in sections 14.41 and 14. 5^ the vertex contribution 
and the self-energy contribution have the same temperature dependence for M2 ^ M\. Moreover, 
as discussed in sections 14.41 and 14.51 the asymmetries are the same when we exchange the internal 
and the external lepton, therefore the asymmetry for a plus-mode external lepton combined with a 
minus-mode internal lepton is the same as the asymmetry for a minus-mode external lepton with a 

plus-mode internal lepton, in short, A7^ = A7 y. We see that generally, the thresholds are the 

ones we expect from our analysis of the decay rates in section [3l For the one-mode calculations we 
have the expected thresholds at for mi and at for \f2mi. For all asymmetries where a plus- 
mode lepton is involved, that is A7_|__|_, A7^ and A7 |-, the phase space is reduced similar to the 

\/^m£ case below T_j^ and an additional reduction of the phase space sets in between and 
since large momenta k or k' that correspond to a large mass m{k) become kinematically forbidden. 
Between these two thresholds, and , the asymmetry for A7_| / becomes larger than the 
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asymmetry for A7_|__|_. This effect occurs because in the (++)-asymmetry the phase spaces of both 
the internal and the external lepton are suppressed, while for the mixed modes, (H — ) or ( — h), only 
the phase space of one momentum is suppressed, while the phase space of the other momentum 
is still large. The effect is similar to the observation that 7_ becomes larger than 7+ above . 
Relying solely on phase-space arguments, one would expect that ^7^^^ is a good approximation 
for A7++. The fact that A7++ is clearly smaller than A7^^^ is due to two suppressing factors: 
One factor is the effect of the two residues Zh{k) and Zh'{k'), which suppress the rate somewhat for 
small momenta k and k' . The other, more important factor is the fact that the helicity structure 
and angular dependence of the integrals are different for the (++)- and the y/2 m£-case as explained 
in section 14. 6i Since neutrino momenta are of the order ~ Mi ~ T for our temperature range, the 
lepton momenta will be of the same order, that is A; > rrii, and the leptons and Higgs bosons will 
preferentially be scattered forward. Thus also the angle ^ between the two leptons will be small 
and the factor H_ defined in equation (jl06p is suppressed, while the corresponding one-mode factor 
1 — ^{kk')/{uJkUJk') is larger than H- for small angles and still finite if both leptons are scattered 
strictly in the same direction, that is = 1. We have checked numerically that this is the main 
reason why A7^^^ > A7_|_+ in the range 1/2 Mi ^ T < . 

Since the CP-asymmetries follow the corresponding finite-temperature decay rates that are 
shown in figure it is very instructive to normalise them via these decay rates, that is 7+, 7_, 7^^ 
and 7^^^- This also gives a more intuitive definition of the CP-asymmetries at finite temperature. 
These asymmetries are shown in figure [T2l normalised by the zero temperature CP-asymmetry. 
Compared to the normalisation via 70 in figure \TT\ we see that the (-|— |-)-asymmetry does not fall 
as steeply as the corresponding decay rate 7+ between and Tg^, so the ratio A7++/7_|_ is dented 
at the threshold Tj^. This illustrates that the (-|— |-)-CP-asymmetry shows a stronger suppression 
below the threshold , since it suffers from two phase space reductions and two residues that 
are smaller than one. Therefore, the (-|--|-)-asymmetry is not affected as strongly as 7+ by the 
additional suppression above when large momenta k and k' are forbidden and the transition 
over this threshold is smoother than for the decay rate 7+. So 7_|_ falls more steeply than A7+_|_ 
above the threshold and the ratio of the two rates has a dent at . For the )-asymmetry, 
this effect is even stronger, since it is less suppressed than the (-|— |-)-asymmetry above , so the 
ratio A7_| /(eo7-) rises up to a value of 0(0.1). 

The CP-asymmetries in Higgs boson decays at high temperature are shown in figure [T3l nor- 
malised to eo7*°* and we have assumed that M2 ^ Mi,T. The behaviour is similar to the neutrino 

decays, where the ( )-asymmetry is strongly suppressed, while the (H — )-, the ( — h)- and the 

(-l--l-)-asymmetries have a strict threshold at Tq and are suppressed due to the reduced phase 
space between Tq and T^, as expected. The (H — )- and ( — h)- asymmetries are the same and are 
somewhat less suppressed than the (-|— |-)-asymmetry between the thresholds Tq and T^. In our 
approximation in section 14.61 we see that the difference of the matrix elements A|A4p rises as 
T^, so A7 rises as T^, which can be seen in the plot for all finite-temperature asymmetries. The 
one-mode asymmetries A7mf, ^l^me (-l--l-)-asymmetry are very close to each other at 

high temperature. 

We also normalise the asymmetries to the corresponding decay rates in figure O and find that 
they all approach a constant value at high temperature, as is expected since the decay rates also 

rise as T^. The dents in the ratios A7_|_+/(eo7+) and A7_| /(eo7+) with an external plus-mode 

lepton are similar to the ones for the neutrino decays and due to the very strong suppression of 
7+ below the threshold T^. The numerically dominant asymmetries A7m^, A7^^^ and A74_-|_ all 
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Figure 12: The CP-asymmetries in neutrino decays normalised by the CP-asymmetry in vac- 
uum and the corresponding total decay density at finite temperature, that is A7_|_+/ (7^*60), 
A7-+/(7reo), A7+_/(7i°teo), ^l-/{l'°'e^), A^^J{^'-%) and A7^„/(7^„/o), where the 
CP asymmetries A7 are explained in figure dH We choose Mi = 10^° GeV and M2 » Mi. 



settle at a rather high CP-asymmetry, two orders of magnitude higher than at zero temperature, as 
we expect from our approximate calculation for a Higgs boson at rest in section [46l This is partly 
due to a suppression of |A^p which rises as = g^T'^, but mainly due to the larger difference in 
matrix elements A|7Wp for Higgs boson decays. 

5 Conclusions 

We have performed an extensive analysis of the effects of HTL corrections on CP-asymmetries in 
neutrino and Higgs boson decays. This implies capturing the effects of thermal masses, modified 
dispersion relations and modified helicity structures. We put special emphasis on the influence of 
the two fermionic quasiparticles, which show a different behaviour than particles in vacuum, notably 
through their dispersion relations, but also the helicity structure of their interactions. Our work is 
thus similar to the work done in reference [H] , where the authors of the latter work did not include 
the effects of fermionic quasiparticles and get a different result for the CP-asymmetries, which are 
crucial for the evolution of the lepton asymmetry. We also approximate the lepton propagators by 
zero temperature propagators with the zero temperature mass replaced by the thermal lepton mass 
or the asymptotic mass. We refer to these cases as one-mode approach and compare the thermal 
cases to the zero-temperature case. 

We have calculated HTL corrections to neutrino decays in reference [23j and shortly presented 
the Higgs decay rate at high temperature in reference |44j . while we analyse this decay in detail in 
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Figure 13: The CP-asymmetries in Higgs boson decays normalised by the CP-asymmetry in vacuum 
and the total decay density in vacuum, A7/(7Q°*eo), where the asymmetries A7 are explained in 
figure El We choose Mi = 10^° GeV and M2 > Mi. 
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Figure 14: The CP-asymmetries in Higgs boson decays normalised by the CP-asymmetry in vac- 
uum and the corresponding total decay density at finite temperature, that is A7++/(7^?*eo), 
A7-+/(7^*eo), A7+-/(7l°*eo), A7__/(7l°*eo), A7„^,/(7:„1eo) and A7^„/(7^„/o), where the 
CP asymmetries A7 are explained in figure [TTl We choose Mi = 10^^ GeV and M2 ^ Mi. 
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section [3l While the zero-temperature fermion mass vanishes in the unbroken electroweak phase, 
the resummation of HTL fermion self-energies results in an effective fermion propagator that does 
not break chiral invariance and is split up in two helicity modes. The external fermion states 
therefore behave conceptually different from the ones with chirality-breaking thermal masses that 
have been inserted in the kinematics by hand. Moreover, one has to take care of one additional 
mode, which has implications for the Boltzmann equations [22] . 

We calculated the CP-asymmetries in section HI To our knowledge, this is the first calculation 
of a CP-asymmetry in leptogenesis that includes HTL-corrections and takes into account corrected 
thermal distributions for the particles in the loopl^ We present rules for the product of spinors 
that are related to the fermionic quasiparticles in equations (|49p and (j47p and derive frequency 
sums for the HTL fermion propagator in equation (I99p . We find four different CP-asymmetries 
corresponding to the four different choices of lepton modes both in the loop and as external states. 
We find the CP-asymmetry to be symmetric under an exchange of the lepton mode in the loop 

and the external lepton mode, such that = A7 ^. At finite temperature, there are three 

possible cuttings for the vertex contribution, the {£' ,<j)'}-cut that corresponds to zero temperature 
and two additional cuts involving the internal namely through and {N',(p'}, which 

have been found by references |14lll91 [l5] and examined more closely in reference [20], using the 
real-time formalism. We obtain the same cuts using the imaginary time formalism and concentrate 
on the {£' ,(p'}-cut, assuming the hierarchical limit of M2 ^ Mi. As expected from the zero- 
temperature result, we find the vertex contribution proportional to the self-energy contribution in 
this limit. Contrary to reference [13], we find that the CP-asymmetry in Higgs boson decays is 
larger than the asymmetry in neutrino decays by about a factor of 100 if appropriately normalised. 
This is due to a suppression of the Higgs boson decay rate and a thermal enhancement of the 
CP-asymmetry by the distribution functions of the Higgs bosons and leptons. We compare the 
CP-asymmetries in the two-mode approach to the CP-asymmetries in the one-mode approach. We 
find that for the two-mode approach, the helicity structure of the modes prohibits the two leptons 
to be scattered strictly in the same direction while for the one-mode approach, this direction is 
only mildly suppressed. Notably this fact is responsible for suppressing the (+-|-)-CP-asymmetry 
compared to the asymmetries of the one-mode approach, as well as the residues of the plus-modes 
to less extent. 

Summarising, we argue that for an accurate description of medium effects on leptogenesis, 
the influence of thermal quasiparticles, notably the effects of the two fermionic modes, cannot be 
neglected. The corresponding CP-asymmetries show a considerable deviation from the one-mode 
asymmetries in the interesting temperature regime T ~ Mi. Moreover, the presence of a new 
minus-mode that essentially does not interact with the SM turns out to have implications for the 
lepton asymmetry even in the strong washout regime, as we found in an additional study |22j . 

^'^Reference [H] has a calculation for the CP-asymmetry of the neutrino and Higgs decays with HTL corrections 
in the propagators and thermal masses in the external states. They obtain a different combination of thermal 
distributions for the particles in the loop. The discrepancy is discussed in detail in reference [21]. The authors of 
the first reference get a factor 1 — fii + f^' — ft ftp' for the neutrino decays and /^/ — ~ S/^'/f for the Higgs 
boson decays due to an erroneous choice of cutting rules in the real time formalism. The correct calculation gives 
1 ~ /f + /(p' for the neutrino decays and /^z + f^i for the Higgs boson decays. This discrepancy is also responsible 
for our CP-asymmetry in Higgs decays being a factor ten larger than their result. 



27 



Acknowledgements 



We would like to thank Mathias Garny, Georg Raffclt, Michael A. Schmidt and Markus Thoma 
for their support and comments in this project. Thanks also to Denis Besak, Dietrich Bodeker, 
Wilfried Buchmiiller, Valerie Domcke, Marco Drewes, Andreas Hohenegger, Alexander Kartavtsev 
and Christoph Weniger for fruitful and inspiring discussions. 



A Frequency Sums 

A.l Frequency sums for HTL fermion propagators 

In order to deal with the HTL lepton propagator, we derive frequency sums for the propagator 
parts A±{K) of a fermion propagator. We write the propagator in the Saclay representation as 

AhiK) = - rdre'=°-A^T,k), 
Jo 

POO 

Ah{T,k)= du phfFi-uj)e-'^\ (94) 



where fp stands for a Fermi-Dirac distribution. Since we are only interested in the pole contribution, 
we write the corresponding spectral density as 

^pole ^ -Zf^[5{uj - COh) + S{UJ + OO^h)] = - X] ^shS{(^ - SUJsh) , (95) 

s 

where sh in Zgh and 0sh denotes the product of s and h, that is, Zgh = for s = /i = — 1 for 
example. We have for the propagator 

/oo 
da; X Zshd{u; - su:sh)fF{-co)e-^^ 
-oo g 

= - E ZshfF{-sUsH)e-^-^'^^ = J2 

s s 

(^) = E ZshfFi-sUs,) f dr eC^o— = X Alf{K) , 

s s 



■'sh 



ko - SLOsh ' 



(96) 



where Zgh = (w^/j — ^^)/(2?Ti^) is the quasiparticle residuum. 

In dealing with frequency sums of bare thermal propagators, it is very convenient to write 

As{K) = A_s{-K). (97) 

Replacing a boson by a fermion amounts to replacing fsi^) by — /f('^)- Moreover, calculating a 
frequency sum of ko times the propagators amounts to replacing ko with so; as in 

T X koAs, {K)As, {P-K) = siooT ^ A,, {K)As, {P - K), (98) 

ko ko 
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where oj = ^k? + m? and m is the mass of the first boson. The same holds for fermions. 

It is straightforward to work out the frequency sums for the resummed lepton propagator, 

^ ' ' 2w po - SlOJs^h - S20J 

feo 

where the other necessary frequency sums can be derived from this by making the appropriate 
substitutions. 

A. 2 The frequency sum for the vertex contribution 



The frequency sums for the pole part of a HTL fermion propagator is derived in appendix lA.il We 
calculate the frequency sum of the three propagators in the vertex loop by partial fractioning 



We are using /\n'[P') = ^n'{-P') and 
The frequency sum is given by 



-a 



^shN' 



sh4>' 



N. 



shN' 



where 



^shN 



/ = 1 - fpisUJsh') - frisN'^^N') , 



and 



^shN' = QO - SUsh' - Sn'UJn' , 
Nsh<l>' = Po — SUJsh' — S^'LO^i . 

Summing over all propagator parts and the helicity-over-chirality ratios, we get 



k'o h' h' P 



{E^H^ + E+H^} , 



where 

H± = l± hh'^ , 
^ = k • k' . 

The coefficients E± are given by 

E^ = F^'^Af - F^'^'Af - F^^'A" + F^'^'Af, 
E+ = F^''^'Af - F^^'4 - F^'^'A^, + F^'Af, 



(100) 



(101) 



(102) 



(103) 



(104) 



(105) 



(106) 

(107) 
(108) 
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and the coefficients F^^ read 



p4,N 


= Bt 




pi'N 


- R^' 


- Bn^ 




= Bi 


- Bn^ 




-Bs 


R^' 
- Bn^ 



?N'0 



S5 



AT' 



rO 
rO 

Bn- 



The factors -Bat/,/, and are given by 



7^ 



N 



N/4, 



N 



where the numerators and denominators read 



AT; 



AT 



AT 



N 



N 



N] 



N 



N' 



N 



Pq + u;' + u;q>, 
Pq-u;' + u;q>, 

Pq+u;' -Uq', 



AT? 



N, 



N' 



ko — ^q' — ^p' , 
fco + LOqi + Wp/, 
fco — Ulqi + Wp/ , 



= qo + ^' + ^p', 



Ni 



and 



-(1-/^' + /</>'), 



<p 



= qo + i^' - 

! — — /at' , 
-(1 - fi' — /at') , 
-{fc - In'), 
fe' — In' ■ 



We can write 



Zh' 



+ 



or, more explicitly, 



(109) 



(110) 



(111) 



(112) 



(113) 



+ [(ijf - <) - - 4 - (i?;' - A" + (B° - B%) Af] H^} . 

(114) 
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B Analytic Expressions for the CP- Asymmetries 
B.l Vertex cut through {£', </>'} 

We simplify the analytic expression for e^(T) in equation ([35]) . For ly in equation (f53|) we get 

= ^^Y. r^k' rdc^'-ZH^Z^iAi-Af)H., (115) 

where it is sufficient to integrate (j)' from to vr. The difference of the matrix elements reads for 
the vertex contribution 



2] MiM2ZhOJh 



4^2 



\M{N ^ 40)1' - \M{N ^ 40)1' = - 5s(/(2)Im | [(a+a) 

X V / dk' d^' — Zh'Z^{A[ - Af )H_ . (116) 
^Jo Jo ^P' 

Correspondingly, the difference in decay rates reads 
7(iV ^ 40) - liN ^ 40) = - 55C/(2)Im I [(At a) 



2] M1M2 



X 



4(27r)5 

V / dEdkdk' r d4>'kFNhZh—Z^Zh'{4 - Af)H^ , (117) 
^ J Jo P0p' 



where F]\rh = /^(l + /^^)(1 — fih) statistical factor for the decay. 

We know from section [3] that 

^ \MUN ^ LH)\^ = gsui2)9c{\^^)iiZhUj{P0 - hpv) , (118) 

s 

where gc = indicates that we sum over N ^ (pi and A'" — )• 0£. Thus 

T{N ^ LhH) = 55^(2)5e-^^^ j dkkZoZhipo - hpv) (119) 

and 

7(A^ ^ LhH) = 55^(2)5, I dEdkkfNZoZhiPo - hpr,) , (120) 

where we have summed over the neutrino degrees of freedom. 
We arrive at 

eH{T) = - 9sui2) ^ L j d^^dMfc d0 kF^J,—Z^Z„{A, - A, )if_ 

Im{[(AtA)i2]n / dii^dMfc- /; d4>'kF^^J,^Z^Zh.{A[ - Af)H^ 

9,(AtA)ii 47r2 JdEdkkfMZoZhipo-hpr]) ' ^^^^^ 
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B.2 Self-energy cut 

For the self-energy contribution, we get 



The difference in decay rates reads 

M1M2 



X 

h- 



{Ml - M|) (27r)5 
I AEdkdk' I d^'kF'^^Zh-Z^Zh'H^ . (123) 



The CP-asymmetry reads 



' -^{N ^ LhN) Ml - Ml{2-Kf ^ J Jo Nu p N h 

Im{[(AtA)i2]n ^1^2 1 T.h' ! ^Edkdk' j; d<P'kF''^lZt,^Z^ZyH. 

<7c(AtA)ii M2-M|27r2 J dEdkkfNZDZh{pQ - hpri) ' ^ ' 

C The Other Cuts 

C.l Imaginary Parts 
Vertex cut through {N2,cf)'} 

We use the conventions for the vertex contribution in A^-decays in section HI For A'^' , we shift 
integration variables to d'^q' after carrying out the Matsubara sum over k'^. We consider the angle 

W = ^ (125) 



between k and q' and write 



where the angle is 



and 



flkq'fi = {-2u!U!q> + Sfc) , (127) 



Sfe = M| - (^2 _ _ ^2 _ ^^28) 
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The imaginary part reads 



— ^A^,A^,A,,i/_) = -31m T / dgV^dv/ d</>,,A^,A^,A,,i/. 



Vertex cut through {N2,£'} 

For the A^^'-term, we integrate over k' , and choose the polar angle between q and k', 

qk' 



(129) 



Vqk' 



qk' 



We write 



where the angle is 



and 



The imaginary part is given by 



-vr 



qk 



Im T 



rim T 



fen 1^' 



(130) 



(131) 



(132) 



(133) 



/•TT 

Jo 



Y^EJ dk'd<p'^Z,,Z^'iAf-A'^,)H^. (134) 



C.2 Analytic Expressions for the CP- Asymmetries 
Vertex cut through {N2,(p'} 

We get for Nf 

M1M2 ZhUJ 



/"OO /"TT / 

16.2 * i: X 1 d0,,A-Z„ [{b'; - Bi,) H. + (b« - <) ff,] . (135) 
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The difference in decay rates reads 

7(iV ^ ih4>) - l{N Ih4>) = -55C/(2)Im 



AtA 



12. 



2] M1M2 



hi •' 



h' 



'0 "<P."^-iv,-^^-^' 



4(27r)5 



+ 



(136) 



and the CP-asymmetry reads 

Im{[(AtA)i2]2}MiM2 

e/ii-' j = - gsu(2) 



I dEdkdq' £ d<P,,kF'^^Zhj^Zh, 



7(iV ^ L/jTV) 4(27r)5 
(5j - B^) + (so - Bt:) K 

Im{[(AU)i2]-}MiM2 
5c(AtA)ii 47r2 

EH'IdEdkdk'J^dcl>'kF^lZf,J^Zf,,[{Bi; - B%)H^ + (5° - 



J dEdkkfNZoZhipo - hpij) 



(137) 



Vertex cut through {Ar2,^'} 

For N^' , we get 



Im(-^v)jv^' 
The difference in decay rates reads 

7(A^ ^ 4</>) - 7(-^ ^ 4^) = - 55C/(2)Im 



AtA 



12. 



2] M1M2 
4(27r)5 



V / dEdkdk' r d<l>gk'kF'^lZh^zf Zh'iAf - A''i)H+. 

Jo Q^q' 



(139) 



The CP-asymmetry reads 



Im{[(At A)i2]n M1M2 E.' / dEdkdk' /; d0,,,A;P^^Z,^zf Z,,(< - ^?)i/- 



5e(AtA)ii 47r2 



J dEdkkfNZoZhipo - hpij) 



(140) 



where we integrate over (f)^^' and we take the coordinate system differently than for N^. 
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D Approximation for the One-Mode Approach at High Temper- 
ature 



Let us examine the high temperature behaviour of the one-mode approach by calculating the CP- 
asymmetry in the matrix elements of a Higgs boson at rest, where we assume that Mi,m£ <C 
nifj, <C M2. For simplicity, we calculate the self-energy contribution. The integral that corresponds 
to equation (|78]) reads 

j^M2Mi [Aat, A^, A^,] *K.K'. (141) 

The part that contributes to the imaginary part of the diagram is 

lo/ilti = 2 / -|^M2MiA^,-^<(K • K') , 

Im< = -7rZ^'5(<) = -vrg;Z^'(5(e - Co) , (142) 



where ^ = {kk') /{kk'), 



^0 = (143) 



and we have neglected Mi and m^. We get 

lm{IorX = -^I^ dk'zt;i2k' - m^) , (144) 
where x = M|/M^ and k = For simplicity, we make the approximation 

4' = A' + h « e-"^"' + e--^'^ = (1 + , (145) 

so 

iMhlDi = -^T^(l + e"'^) r ^k'e-^'^{2k' - m^) . (146) 
Svr 1 — X jfc/ 

The integral gives 

Im(/o/r)t: = -l^r^(l + e-''f^)T\-'^ . (147) 
47r 1 — X 

We parameterise as = g^T and obtain 



A = -\M\'^ = -8ImAcpIm(/o/n = 2 ^™^^^^^ rV^^/^(l + e'^^/^^ . (148) 

vr 1 — X 

Using the expression 

\Mtot\^ = 4(AtA)nK • P = 2(AtA)ii5^r2 (149) 
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we arrive at 



IMotr ('^^A)ii7rl-X5^ 



Assuming that (7^ ^ 1, we get 



- ^ (151) 



Taking f/,^ = m^jy/T ~ 0.42 for T = 10^^ GeV and using the more accurate term in equation (|150p . 
we get e/eo ~ 70, while we get e/eo ~ 90 for equation ()15ip . We view this result as a rough 
approximation of the value of the CP-asymmetry in Higgs boson decays at high temperature. Both 
our approximation and the numerical solution of the exact expression in section 14.71 give a factor 
of 100 difference to the CP-asymmetry in vacuum. 



References 

[1] E. Komatsu, K. Smith, J. Dunkley, C. Bennett, B. Gold et. al, Seven-Year Wilkinson 
Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation, 
[arXiv: 1001.4538]. 

[2] M. Fukugita and T. Yanagida, Baryogenesis Without Grand Unification, Phys. Lett. B174 
(1986) 45. 

[3] P. Minkowski, — )• 67 at a Rate of One Out of 1-Billion Muon Decays?, Phys. Lett. B67 
(1977) 421. 

[4] T. Yanagida, Horizontal Gauge Symmetry and Masses of Neutrinos, In Proceedings of the 
Workshop on the Baryon Number of the Universe and Unified Theories, Tsukuba, Japan, 
13-14 Feb 1979. 

[5] M. Gell-Mann, P. Ramond and R. Slansky, Complex Spinors and Unified Theories, 
Print-80-0576 (CERN). 

[6] R. N. Mohapatra and G. Senjanovic, Neutrino Masses and Mixings in Gauge Models with 
Spontaneous Parity Violation, Phys. Rev. D23 (1981) 165. 

[7] J. Schechter and J. W. F. Valle, Neutrino Masses in SU(2) x U(l) Theories, Phys. Rev. D22 
(1980) 2227. 

[8] J. Schechter and J. Valle, Neutrino Decay and Spontaneous Violation of Lepton Number, 
Phys.Rev. D25 (1982) 774. 

[9] F. R. Klinkhamer and N. S. Manton, A Saddle Point Solution in the Weinberg- Salam 
Theory, Phys. Rev. D30 (1984) 2212. 

[10] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, On the Anomalous Electroweak 
Baryon Number Nonconservation in the Early Universe, Phys. Lett. B155 (1985) 36. 



36 



[11] A. D. Sakharov, Violation of CP Invariance, C Asymmetry, and Baryon Asymmetry of the 
Universe, Pisma Zh. Eksp. Teor. Fiz. 5 (1967) 32-35. 

[12] S. Davidson, E. Nardi and Y. Nir, Leptogenesis, Phys. Kept. 466 (2008) 105-177 
|arXiv:0802.2962j. 

[13] L. Covi, N. Rius, E. Roulet and F. Vissani, Finite Temperature Effects on CP Violating 
Asymmetries, Phys. Rev. D57 (1998) 93-99 |hep-ph/9704366|. 

[14] G. F. Giudice, A. Notari, M. Raidal, A. Riotto and A. Strumia, Towards a Complete Theory 
of Thermal Leptogenesis in the SM and MSSM, Nucl. Phys. B685 (2004) 89-149 
|hep-ph/0310123| . 

[15] A. Anisimov, W. Buchmiiner, M. Drewes and S. Mendizabal, Quantum Leptogenesis I, 
|arXiv: 1012.5821] . 

[16] A. Anisimov, D. Besak and D. Bodeker, Thermal Production of Relativistic Majorana 
Neutrinos: Strong Enhancement by Multiple Soft Scattering, |arXiv: 1012. 37M| . 

[17] M. Garny, A. Hohenegger and A. Kartavtsev, Quantum Corrections to Leptogenesis from the 
Gradient Expansion, [ arXiv: 1005.5385] . 

[18] M. Beneke, B. Garbrecht, C. Fidler, M. Herranen and P. Schwaller, Flavoured Leptogenesis 
in the CTP Formalism, Nucl. Phys. B843 (2011) 177-212 [arXiv: 1007.4783] . 

[19] M. Beneke, B. Garbrecht, M. Herranen and P. Schwaher, Finite Number Density Corrections 
to Leptogenesis, Nucl. Phys. B838 (2010) 1-27 |arXiv: 1002. 1326] . 

[20] B. Garbrecht, Leptogenesis: The Other Cuts, [arXiv: 1011 .3122]. 

[21] M. Garny, A. Hohenegger and A. Kartavtsev, Medium Corrections to the CP-Violating 
Parameter in Leptogenesis, Phys. Rev. D81 (2010) 085028 [arXiv: 1002 . 0331] . 

[22] C. P. Kiefiig and M. Pliimacher, Hard- Thermal- Loop Corrections in Leptogenesis U: Solving 
the Boltzmann Equations, [arXiv: 1111 . 1235] . 

[23] C. P. Kiefiig, M. Pliimacher and M. H. Thoma, Decay of a Yukawa Fermion at Finite 
Temperature and Applications to Leptogenesis, Phys. Rev. D82 (2010) 036007 
[arXiv: 1003.3016] . 

[24] M. Le Bellac, Thermal Field Theory. Cambridge University Press, Cambridge, 1996. 

[25] E. Braaten and M. H. Thoma, Energy Loss of a Heavy Fermion in a Hot Plasma, Phys. Rev. 
D44 (1991) 1298-1310. 

[26] E. Braaten and M. H. Thoma, Energy Loss of a Heavy Quark in the Quark- Gluon Plasma, 
Phys. Rev. D44 (1991) 2625-2630. 

[27] O. K. Kalashnikov and V. V. Klimov, Polarization Tensor in QCD for Finite Temperature 
and Density, Sov. J. Nucl. Phys. 31 (1980) 699. 



37 



[28] D. J. Gross, R. D. Pisarski and L. G. Yaffe, QCD and Instantons at Finite Temperature, 
Rev. Mod. Phys. 53 (1981) 43. 

[29] E. Braaten and R. D. Pisarski, Soft Amplitudes in Hot Gauge Theories: A General Analysis, 
Nucl. Phys. B337 (1990) 569. 

[30] E. Braaten and R. D. Pisarski, Deducing Hard Thermal Loops From Ward Identities, Nucl. 
Phys. B339 (1990) 310-324. 

[31] H. A. Weldon, Effective Fermion Masses of Order gT in High Temperature Gauge Theories 
with Exact Ghiral Invariance, Phys. Rev. D26 (1982) 2789. 

[32] E. Braaten, R. D. Pisarski and T.-C. Yuan, Production of Soft Dileptons in the Quark - 
Gluon Plasma, Phys. Rev. Lett. 64 (1990) 2242. 

[33] E. Braaten and R. D. Pisarski, Galculation of the Quark Damping Rate in Hot QCD, Phys. 
Rev. D46 (1992) 1829-1834. 

[34] V. V. Klimov, Spectrum of Elementary Fermi Excitations in Quark Gluon Plasma. (In 
Russian), Sov. J. Nucl. Phys. 33 (1981) 934-935. 

[35] R. D. Pisarski, Renormalized Gauge Propagator in Hot Gauge Theories, Physica A158 
(1989) 146-157. 

[36] J. I. Kapusta, P. Lichard and D. Seibert, High-Energy Photons from Quark - Gluon Plasma 
Versus Hot Hadronic Gas, Phys. Rev. D44 (1991) 2774-2788. 

[37] F. Karsch, M. G. Mustafa and M. H. Thoma, Finite Temperature Meson Correlation 
Functions in HTL Approximation, Phys. Lett. B497 (2001) 249-258 [hep-ph/0007093|. 

[38] J. M. Cline, K. Kainulainen and K. A. Olive, Protecting the Primordial Baryon Asymmetry 
from Erasure by Sphalerons, Phys. Rev. D49 (1994) 6394-6409 |hep-ph/9401208|. 

[39] P. Elmfors, K. Enqvist and I. Vilja, Thermalization of the Higgs Field at the Electroweak 
Phase Transition, Nucl.Phys. B412 (1994) 459-478 [he p-ph/93072 10]. 

[40] K. Kajantie, M. Laine, K. Rummukainen and M. E. Shaposhnikov, Generic Rules for High 
Temperature Dimensional Reduction and their Application to the Standard Model, Nucl. 
Phys. B458 (1996) 90-136 |hep-ph/9508379| . 

[41] H. A. Weldon, Simple Rules for Discontinuities in Finite Temperature Field Theory, Phys. 
Rev. D28 (1983) 2007. 

[42] R. L. Kobes and G. W. Semenoff, Discontinuities of Green Functions in Field Theory at 
Finite Temperature and Density. 2, Nucl. Phys. B272 (1986) 329-364. 

[43] C. P. Kiefiig and M. Pliimacher, Thermal Masses in Leptogenesis, AIP Conf. Proc. 1200 
(2010) 999-1002 [arXiv: 0910 .4872) . 

[44] C. P. KieBig, M. Pliimacher and M. H. Thoma, Fermionic Quasiparticles in Higgs Boson and 
Heavy Neutrino Decay in Leptogenesis, J.Phys.Conf.Ser. 259 (2010) 012079. 



38 



[45] M. Garny, A. Hohenegger, A. Kartavtsev and M. Lindner, Systematic Approach to 

Leptogenesis in Nonequilibrium QFT: Vertex Contribution to the CP- Violating Parameter, 
Phys. Rev. D80 (2009) 125027 [arXiv: 0909. 1559] . 



39 



